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Abstract. The inclusion of nonadiabatic corrections to the electron-phonon interaction leads to a strong 
momentum dependence in the generalized Eliashberg equations beyond Migdal's limit. For a s-wave sym- 
metry of the order parameter, this induced momentum dependence leads to an enhancement of T c when 
small momentum transfer is dominant. Here we study how the rf-wave symmetry affects the above behavior. 
We find that the nonadiabatic corrections depend only weakly on the symmetry of the order parameter 
provided that only small momentum scatterings are allowed for the electron-phonon interaction. In this 
situation, We show that also for a d-wave symmetry of the order parameter, the nonadiabatic corrections 
enhance T c . We also discuss the possible interplay and crossover between s- and d-wave depending on the 
material's parameters. 

PACS. 63.20.Kr Phonon-electron and phonon-phonon interactions - 71.38.+i Polarons and electron- 
phonon interactions - 74.20.Mn Superconductivity: nonconventional mechanisms 



1 Introduction 

In ordinary low-temperature superconductors, the small- 
ness of the relevant phonon frequency <x>o compared to the 
Fermi energy Ep permits to formulate a theory of super- 
conductivity based on a closed set of formulas known as 



allowing quanti- 
The closed form 



Migdal-Eliashberg (ME) equations ] 
tative agreements with experiments 
of the ME equations stems from Migdal's theorem which 
states that as long as ivo/Ep -c 1 the electron-phonon (e- 
ph) vertex corrections to the electron self-energy are at 
least of order \u>o/Ep, where A is the e-ph coupling, and 
can therefore be neglected 

A different situation is encountered in high-T c super- 
conductors such as cuprates and fullerides. These materi- 
als have in fact Fermi energies much smaller than those 
of conventional metals j|j5] so that the energy scale ujq 
associated to the mediator of the superconducting pairing 
can be comparable to E-p. Hence, the quantity luq/Ef is 
no longer negligible and in principle vertex corrections be- 
come relevant preventing the ordinary ME scheme to be 
a correct description of the superconducting state. 

The possible breakdown of Migdal's theorem in high- 
T c superconductors inevitably calls for a generalization 
beyond the ME scheme to include the no longer negligible 
vertex corrections. A possible way to accomplish this goal 
is to rely on a perturbative scheme by truncating the infi- 
nite set of vertex corrections at a given order. In previous 
works, we have proposed a perturbative scheme in which 
the role of small parameter is played roughly by Xluq/Ep 



leading to a generalized ME theory which includes the 
first nonadiabatic vertex corrections [HQ]. For a single- 
electron Holstein model, such a first order perturbative 
approach leads to good agreements with exact results as 
long as the system is away from polaron formation, that 
is for A < A c ~ 1 ||. The region of validity of the per- 
turbative approach, which we could name nonadiabatic 
region, is characterized by quasi- free electrons (A < 1) 
coupled in a nonadiabatic way (coo/Ep not negligible) to 
the lattice. According to this definition, our nonadiabatic 
region is different from the classic polaronic picture. Of 
course, larger values of A would render higher order ver- 
tex corrections important leading to the breakdown of our 
truncation scheme. 

The key point of the nonadiabatic theory, is that the e- 
ph effective interaction is described in terms of vertex cor- 
rections which depend on the momentum transfer |q| = q 
and the Matsubara exchanged frequency lu in a non-trivial 
way. For example, the e-ph vertex correction appearing in 
the normal state self-energy becomes positive (negative) 
when vpq < uj (vpq > to), where vp is the Fermi velocity 
jffl. The generalization to the superconducting transition 
reveals that this situation is also encountered for the class 
of diagrams beyond Migdal's limit relevant for the Cooper 
channel. Concerning the critical temperature T Cl as long 
as the momentum transfer is less than luq/vf, the nonadi- 
abatic corrections lead to a strong enhancement of T c also 
for moderate values of the e-ph coupling A Such 
a strong momentum-frequency dependence of the vertex 
corrections is confirmed by numerical calculations within 
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a tight-binding approach and general theoretical con- 
siderations on the phy sical interpretation of such nonadi- 
abatic corrections H]. 

So far, nonadiabatic superconductivity has been stud- 
ied by requiring the order parameter A to be independent 
of the momenta. This situation is certainly suitable for 
the fullerene compounds which are s-wave superconduc- 
tors. However, one striking characteristic of several high- 
ly superconductors is the strong momentum dependence 
of the order parameter A(k). Among the several types of 
measurements aimed to resolve the pairing symmetry, the 
Josephson tunneling |l2| and angle resolved photoemis- 
sion experiments |13| are the most convincing ones show- 
ing that the order parameter of several cuprates, maybe 
with the exception of the electron-doped NCCO, has a pre- 
dominant d- wave symmetry: A(k) ~ A[cos(k x ) — cos(k y )]. 

The origin of the d-wave symmetry in high-T c cuprates 
is still debated. On one hand, the observed d-wave symme- 
try is regarded as an evidence against any purely electron- 
phonon pairing interaction so that the mechanism respon- 
sible for superconductivity should be sought among pair- 
ing mediators of electronic origin (like antiferromagnetic 
fluctuations) with eventually a minor electron-phonon com- 
ponent. On the other hand, several theoretical studies 
have shown that the e-ph interaction could produce, un- 
der some quite general circumstances, a d-wave symmetry 
of the condensate Jl4|,[l5| ■ This could happen when for ex- 
ample charge carriers experience an on-site repulsive in- 
teraction together with a phonon induced attraction for 
large inter-electrons distances. The on-site repulsion in- 
hibits the isotropic s-wave superconducting response lead- 
ing the system to prefer order parameters of higher angu- 
lar momenta. A quite general analysis of the interplay be- 
tween on-site repulsion and neighbour and next-neighbour 
attraction has shown s-wave to d-wave crossover depend- 
ing on the microscopic parameters of a model BCS Hamil- 
tonian Jl6| . 

The purpose of the present paper is to study how the d- 
wave superconducting response resulting from a strongly 
momentum dependent total interaction, is affected by the 
inclusion of nonadiabatic vertex corrections. In particu- 
lar, we intend to clarify whether the complex momentum- 
frequency structure of the nonadiabatic contributions could 
sustain an underlying d-wave symmetry of the order pa- 
rameter. 

In the next section we introduce the model and the 
corresponding ME equations for s- and d wave symmetries 
of the gap. In Sec. 3 we generalize the ME equations to 
include the nonadiabatic terms for each symmetry chan- 
nel and calculate the corresponding critical temperatures. 
We find that the theory of nonadiabatic superconductivity 
can lead to d-wave even for phonons in a broad parameter 
range which depends on the degree of electronic correla- 
tion. 



2 The model 

In this section we introduce a simple model interaction 
suitable for our investigation beyond Migdal's limit and 
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Fig. 1. Sketch of the total (el-ph + Coulomb) interaction in 
momentum space. 

capable of providing for s- or d-wave symmetries of the or- 
der parameter. Let us consider the anomalous self-energy 
at the critical temperature 

S s (k) = Wfc - k')G(k')G(-k')Z s (k') , (1) 

k' 

where G(fc') is the fermion dressed propagator: 

G ^ = 1 r 7U-, (2) 

and is the normal self-energy. We use the compact no- 
tation k = (k, Lu n ), k' = (k', ui m ) and J2k> = Em J2k' 
where u n , oj m are fermionic Matsubara frequencies and 
k, k' are electronic momenta (from now on, all momenta 
are two-dimensional vectors lying on the copper-oxygen 
plane) . 

To define the model interaction V pa i r (fc — k') we have 
made use of a number of informations gathered from pre- 
vious studies. First, in order to obtain order parameters 
with higher angular momenta than s-wave, it is sufficient 
to consider a pair interaction made of a repulsive part at 
short distances and an attractive one at higher distances 
(Fig. 1). In momentum space, this interaction corresponds 
to an attractive coupling for small q and a repulsive one 
for large q, where q = k — k' is the momentum transfer. 

Let us now try to interpret this strong momentum 
modulation in terms of e-ph and electron-electron interac- 
tions. In strongly correlated systems, the e-ph interaction 
acquires an important momentum dependence in such a 
way that for large values of the momentum transfer q 
the e-ph interaction is suppressed, whereas for small val- 
ues of q it is enhanced |T^]. A physical picture to justify 
this momentum modulation is the following |l8|] . In many- 
electrons systems a single charge carrier is surrounded 
by its own correlation hole of size £ which can be much 
larger than the lattice parameter a in the strongly corre- 
lated regime. This implies that one electron interacts with 
molecular vibrations of wavelength of order £ or larger, 
leading to an effective upper cut-off q c ~ in the mo- 
menta space. Thus we have a non zero electron-phonon 
interaction when |q| < q c . The cut-off momentum q c can 
also be regarded as a measure of the correlation in the sys- 
tem: aq c <C 1 in strongly correlated systems while aq c ~ I 
in the case of free electrons. 
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From the above considerations, the attractive part at 
small q of our model pairing interaction finds a natural 
interpretation in the e-ph coupling modified by the strong 
electron correlations. We introduce therefore the following 
simple form for the e-ph part of the pairing interaction: 

V(k-k') = \g(k--k')\ 2 D(uj n -Lu m ) 



9(q c - \k-k'\)D(uj r , 



where 



D(oj r , 



K - W m ) 2 + UJ 2 



(3) 



(4) 



is the phonon propagator for which we have adopted a 
simple Einstein spectrum with frequency luq- I n Eq. ([}]), 9 
is the Heaviside step function and the prefactor (irkp/qc) 
has been introduced in order to assure that the momen- 
tum average of |g(k — k')| 2 becomes g 2 for relatively small 
values of the cut-off q c regardless of the particular sym- 
metry of the order parameter. In this way the comparison 
between s- and d-wave solutions, especially in the nonadia- 
batic case treated in the next section, is more transparent. 

Having defined the nature of the attractive part of 
the total pairing interaction V^, a i r (fc — k'), we offer now 
a possible interpretation for the remaining repulsive part 
acting at large q. This repulsion is given by the resid- 
ual e-e interaction and its momentum dependence can be 
obtained, in analogy with the renormalization of the e- 
ph interaction, by using the above picture of correlation 
holes. In this picture, the residual e-e interaction should 
ensure that charge fluctuations with wavelength less than 
£ are unfavourable. This can be modelled by requiring 
that in momentum space the residual interaction is repul- 
sive for |k— k'| > q c and, by using again the theta-function 
for later convenience, we introduce therefore the following 
residual repulsion: 



EZ(k-k') = U 



■nkv, 



Ik-k'l-Sc) 



(5) 



In the above expression, U > and the factor Ttk-p/q c has 
been introduced for the same reason as in Eq.(||). Note 
that, in principle, q c depends on U, however here we shall 
treat two quantities independently on each other by keep- 
ing in mind that small values of q c correspond roughly 
to large values of U. An additional simplification of the 
following calculations is achieved by expressing the off- 
diagonal self-energy (|l|) in terms of a suitable pseudopo- 
tcntial U* rather then U . It is then opportune to formally 
replace equation (JsJ) by 



U*(k-k') = U* 



|k-k'|-2c)- 



UJ 2 



(U! n - LU m ) 2 + UJ 2 

(6) 

where U* represents the dynamically screened Coulomb 
repulsion and the last factor is a cut-off over the Matsub- 
ara frequencies which has been chosen to have the same 
functional form of the phonon propagator for convenience. 



By summarizingthe above results, in the off-diagonal 
self-energy Eg, Eq.(|lJ), the total pairing interaction V pll i T (k~ 
k') is given by: 



Vpair(fc ~ k') = V(k - k') +U*{k- k'), 



(7) 



where V(k — k') and U*(k—k') are given by equation s (||) 
and (^J) , respectively. Finally, the normal state self-energy 
En entering (0) is given by 



(8) 



where the electron-electron interaction has been absorbed 
in a shift of the chemical potential. 

In what follows we assume the Fermi surface to be a 
circle in the momenta space; thus the electronic energy 
depends only on |k|. Moreover, we approximate equations 
(||) and (jfy by keeping |k| = |k'| = fcp so that, for ex- 
ample, k = k-p (cos <p, sine/)). In this way both Ss and En 
depend on the momentum k only via the angle <p. At this 
point it is convenient to transform the momentum inte- 
grations appearing in E§ and En into energy integrations 
as follows: 



E 



2tt 



deN(e) 



(9) 



where N(e) is the density of states for the electrons. We 
make the approximation of constant value for N(e) — Nq 
and finite bandwidth E such that the energy is defined in 
the interval —E/2 < e < E/2. The chemical potential is 
fi = 0, so that we refer to the half-filled situations (Ep = 
E/2). 

On performing the integration over the energy, the 
anomalous self-energy E$ reduces to: 



E s (<p,ui n ) = N wT c J2 J ^ [\g(cos0)\ 2 -U*(cos8)] 

E/2 



xD(Ld n - ui m )- — — arctan 

\U m \Z{(p',U m ) 7T 



\LO m \Z(4>', LO m ) 

(10) 



Z{4>,u n ) = 1 -7V — V / ^\g{cose)\ 2 D(uj n -uj m ) 

E/2 



LU m 2 

x r — arctan 



|w m | it 



\uj m \Z(<j)', uj m ) 



(11) 



where E^{^>,uj n ) = iuj n [l — Z(<j),ui n )] and 8 — (f> — 4>' . The 
wave function renormalization Z((f>,ui n ) actually does not 
depend on <j> and reduces to: 



Z(uj n ) = 1 ~N (g 2 ) ^Ey^ D(ui n - w m ) 



x — arctan 

7T 



E/2 



\Um\Z(w m ) 



(12) 
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where (g 2 )o = Jl d8\g(cos8)\ 2 /(2tt). Let us expand the 
off-diagonal self-energy (fLOl) as follows: 



(13) 



£ s (0,u; n ) = ^ 4°K)1W) . 



i=-c 



where Yi(c/>) = e 11 ^ / y/2n are eigenfunctions of the operator 
L = —id/dcj). By requiring -Eg(c/>, w n ) be real and invariant 
under — > ± 7r (singlet pairing) the above expansion 
reduces to: 



^^ + ^4 2) K)cos(20) + ---, (14) 

V27T V 71" 



where we have singled out the s-wave and d-wave compo- 
nents since in the following we consider only these sym- 
metries. By multiplying both sides of ( [To| ) by Yjt{4>) and 
integrating over c/>, it is straightforward to show that the 
equations for different values of the index I are decoupled 

and that Eg\u> n ) reduces to: 



4 Z) K) = - No [{g 2 )i ~ (lT)i] 7rT c 5>(a 



E s ' {w rn ) 2 
x i — — — — arctan 

\UJ m \Z(UJ m ) 7T 



E/2 



where 



l 

27 
1 

2^ 



\uj m \Z{ui m ) 



d9\g(cos 
dOU* (cos 9)e~ lW 



(15) 

(16) 
(17) 



Finally, by introducing the coupling constants A; = 

N {g 2 )i, A*? = N {U*) h and by setting Z\/ = 4°/^ the 
Eliashbcrg equations assume the following more familiar 
form: 



x — arctan 

7T 



S/2 



|w m |Z(w m ) 



Z(uj n )Ai(uj n ) = - (A; - h*)-kT c ^2 D ( L 



x — arctan 

7T 



J5/2 



|w m |Z(w m ) 



(18) 



(19) 



For i = and 1 = 2, equations ( [1 8|) and ( |19|) are 
Migdal-Eliashberg equations for s-wave and d-wave sym- 
metry channels, respectively. The explicit expressions of 
the constants A/ and fi* follow from the models we adopted 
for the electron-phonon and electron-electron interactions. 
For I = (s-wave) they reduce to: 



An = A 



< %c Hk-k'|)) !=0 = A=P^, (20) 



Mo = M 

* 

= M 



9c 



(0(|k-k'|-<?c)> i=o 



7r arcsin Q c 



(21) 



while for Z = 2 (d-wave): 



A, = A 



(%c-|k-k'|)) i=2 
A(1-2Q2)^T— Qf, 



(22) 



M2 = M 



7rfcp 



<0(|k-k'|-?c)) 



/=2 



(23) 



where A = N g 2 , n* = N Q U* and Q c = q c /2k F . Before 
we generalize the above expressions to include the nona- 
diabatic vertex corrections, it is useful to briefly examine 
qualitatively how the the magnitude of the cut-off param- 
eter Q c affects the gap symmetry. The total interaction 
in the gap equation ( |l9|) is weighted by A; — n\. For the 
s-wave channel it reduces to: 



* /a „ N arcsin Qc tt » 
A - Mo = (A + /z)— — „ 

while for the d-wave case Z = 2 it becomes: 



A 2 - M 2 = (A + - 2Q c 2 )Vr^0l 



(24) 



(25) 



When Q c = 1 (that is g c = 2fcp) the repulsive interaction 
([|) vanishes and the e-p/i coupling (|^) becomes structure- 
less. In this limit we expect the s-wave solution to dom- 
inate over the d-wave one. In fact, from ( pi| ) and (p5|), 
Ao — /ig = A and A2 — /i 2 = 0. By lowering Q c , the to- 
tal interaction acquires a momentum dependence, however 
A2 — ^2 remains negative as long as 1 / V2 < Q c < 1 . In this 
range therefore there is not ci-wave solution. By further 
lowering of the cut-off parameter, the d-wave symmetry 
begins to compete with the s-wave one and for Q c <C 1 
A2 — H2 — A + fx* while Ao — /ig — —k[j.*/2Q c signalling 
that the d-wave symmetry overcomes the s-wave ones. 

In previous studies, we have shown that the nonadi- 
abatic corrections lead to an enhancement of the critical 
temperature for small values of Q c in the s-wave channel. 
However in the present model, small values of Q c lead to a 
solution with d-wave symmetry and the question we face 
in the following sections is whether also in this symmetry 
the nonadiabatic corrections provide for an amplification 
of T r . 



3 Non adiabatic vertex corrections 

In this section we introduce the first corrections arising 
from the breakdown of Migdal's theorem in the equations 
for the normal and anomalous self-energies. By following 
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For the momentum dependence in the electron-phonon 
coupling we can approximate 



l9(k-q)| 2 | 9 (q-k 



/\|2 



| 5 (k-k')|l5(k-q)| 2 . (31) 







? ■ 
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Fig. 2. Self-consistent gap equation including the first correc- 
tions beyond Migdal's theorem 

Ref.0], the first nonadiabatic corrections to the e-ph in- 
teraction affect the normal state self-energy (JsJ) in the fol- 
lowing way: 



This approximation is valid for relatively small values of 
the cut-off q c . Therefore the last term in Eq.(p9|) reduces 
to: 

Y, V(k - q)V(q - k')G(q)G(q - k - k') 
1 

~ V(k - k') V{k - q)G{q)G{q -k-k') . (32) 



At this point it is useful to introduce a compact notation 
for the vertex and cross functions: 

P v (k, k') = \j2 V{ - k - l) G (l) G {l - k + k ') > ( 33 ) 



S N (k) = ^2Vu(k,k')G(k'), 



(26) 



P c (k, k') = \Yl V ( k - 9) G (9)G(g - k-k') 



(34) 



V N (k,k') = V(k-k') 



l+^2V(k-q)G(q-k + k')G(q) 



_(27) 

where the last term in the square bracket of Eq.([27|) de- 
fines the vertex function. The off-diagonal self-energy in 
the nonadiabatic regime is instead modified as follows: 

S s (k) = J2[v s (k,k')-U*(k-k') 



x G(k')G(-k')Ss(k'), 



(28) 



V s {k,k') = V(k-k') 



1+^2 v ( k " l)G(q)G(q -k + k') 



Y J V(k-q)G^q)G{-q + k-k l ) 
Y i V{k-q)V{q-k')G{q)G{q-k 



Thus equations ( p7j ) and ( p9| ) may be written in a simpler 
way as follows: 

V N (k, k') = V(k - k') [1 + \P v (k, k 1 )] , (35) 
V s (k, k') = V{k - k') [1 + 2\P v {k, k') + XP c (k, k')] . 

(36) 

In Eq.(p8|), the momentum dependence of Ssik) is trans- 
formed as in Eq.([l3|) and the interaction term Vs(/c, k') is 
replaced by its angular weighted average: 

(V s (k,k'))i = ^J d6V s (cos6)e- m , (37) 

which in terms of averaged vertex and cross corrections is 
expressed as: 



k% 



lven 



(V s (k, k')) l = (V(k - k')) l + 2X(V(k - k')P v (k, k% 

+ X{V(k-k')P c (k,k'))i. (38) 

(29) The first term in the r.h.s. contains only the e-ph inter- 
where q ee (q, oji) and U*(k - k') is given by Eq.(§). The action and thc P^non propagator and it is simply g: 
second and the third terms within the square brackets in ^Y- 

Eq.(p9|) correspond to the first order vertex corrections, r . 

while the last term corresponds to the cross scattering. (V(k— k'))i = g 2 D(uj n — oj m ) 
These non adiabatic terms are shown in Fig. 2 in which we 
are only include these terms that give a finite contribution _ 



<k 



■D(uj n - u m ), 



k'l)}/ 
(39) 



for T = T C . 

In the vertex corrections there is simple one-phonon i— i 
interaction, while in the cross term the sum refers to the ^ re A ' for 1 = and 1 = 2 is S iven in equations (2C) and 
product of both phonon propagators. This product of phonon O ' respectively The second and third terms of (|38|) cor- 
propagators can be approximated as H: respond instead to the momentum averages of the nona- 

u diabatic corrections 



(Uj n - UJl) 2 + UJ 2 (Ul - LU m ) 2 + UJ 2 



(Wti - LU m ) 2 + UJq (uj n - Ul) 2 + CJq 



(30) 



(V(k - k')Pv(k,k'))i = g 2 D(u> n - u m )Py(uj ni uj mi q c ), 
(V(k ~ k')P c {k, k'))i = g 2 D{cj n - uj m )P^{Lu n ,uj m ,q c ), 

(40) 
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0, u„ 



■ LU 



eters lo„ = 0, w„ 



: LOO 



where 

Py(uj n ,uj m ,q c ) 

P l G (uj n ,u) m ,q c ) 



<7c 

■nk-p 



(% c -|k-k'|)Pv(fc,fc'))i 

(% c -|k-k , |)P C (*,* , ))i. 



(41) 

Analytic expressions of the vertex and cross functions to- 
gether with their averages Py and Pq for I = and 
I = 2 are reported in Appendix and the results are shown 
in Figs. 2 and 3 as function of the adiabatic parameter 
u>o/E-p and for different values of dimensionless cut-off 
Qc = 9c/2fcp. All the curves have been obtained by set- 
ting lu„ — and uj m — ojq so that the exchanged frequency 
equals ojq. The behaviors, particularly at small Q c , of Py 
and Pq are essentially independent of the particular sym- 
metry. In fact for both I = (s-wave) and I = 2 (d-wave) 
the nonadiabatic corrections are positive leading to an en- 
hancement of the total e-ph interaction. We expect there- 
fore that, as for the s-wave case also for the d-wave 
symmetry the vertex and cross corrections tend to amplify 
T c when Q c is sufficiently small. 

To verify this point, we can write down the nonadi- 
abatic Eliashberg equations for different symmetry chan- 
nels. As in the previous section, the normal state self- 
energy ( p6| ) is averaged over the Fermi surface and, ac- 
cording to ([35]), Vn(&, k') reduces to: 

Vn(M') <Vn(M'));=o ■ 

D(uj n - LJ m ) r , 2l ,t0 



No 



[Ao + \ 2 Py-°{uj„,LO m ,qc)] , 



and Z(uj n ) = 1 — S^\uj n ) / iio n becomes: 



Z(w n ) = 1 V [A + A 2 P V =°K, u mi q c j\ 

W n 

m 

E/2 



xD(u n — uj m ) ", — arctan 



\u! m \Z(u} m ) 



(43) 



Finally, the gap function for different symmetry channels 
is: 

Z(uj n )Ai{w n ) = -7tT c ^[A ( + 2A 2 P v (u„,w m ,g c ) 

AW 



+A P G (Lu n ,Lu m ,q c ) - fJ,i]D(iJ n - LO m )- 



x — arctan 

7T 



E/2 



\u> m \Z(u m ) 



(44) 



(42) 



where ^ is given by Eqs.( pi| , p3] ) and I — 0, 2. 

To establish the range of Q c values in which the d-wave 
symmetry is more stable than s-wave one and to quan- 
tify the effect of nonadiabaticity, we solve numerically the 
generalized Eliashberg equations (43) and (|44| ) for I = 
and I = 2. To find T c , we follow the maximum eigenvalue 
method described for example in 0. The resulting values 
of T c as a function of Q c are shown in Fig. 4 for the el- 
and s-wave symmetries. In the inset we show the criti- 
cal temperature calculated without the vertex and cross 
corrections, i. e., for the ME equations ( |l8| ) and (|l9[). To 
display the crossover between the d- and the s-wave sym- 
metries more clearly, we show the results for A = 1 and 
/1* = 0.1. When /i* ~ A in fact the s-wave symmetry is 
suppressed by the strong repulsive interaction. 
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Fig. 5. Behaviour of the critical temperature T c as function of 
Q c in the s- and d-wave symmetry channels. The case shown 
refers to the parameters A = 1 and n* — 0.1. In the inset it is 
shown the case without nonadiabatic corrections. 



By reducing Q c , the s-wave solution (dashed lines) de- 
creases monotonically when the vertex and cross correc- 
tions arc not included (inset) while for the nonadiabatic 
case the corresponding T c shows an upturn before falling 
to zero at Q c — * 0. This latter feature is due to the nonadi- 
abatic corrections which become more positive when Q c is 
small. For lower values of Q c , however, the pseudopoten- 
tial is dominant and T c falls rapidly to zero. Contrary to 
the isotropic case, the d-wave solutions (solid lines) lead to 
critical temperatures which increase when Q c is lowered. 
Since, as discussed before, the vertex corrections have a 
similar behavior both in d- and in s -wave symmetries 
when Q c is small, the critical temperature in the nonadi- 
abatic case is enhanced compared to the solution without 
vertex and cross corrections. 

It is finally interesting to compare the present results 
with the phenomenology of the superconducting copper- 
oxides, which show d-wave, and the fullerides, which in- 
stead show s-wave. In our perspective there are important 
differencies between the two materials. A relevant one is 
that the oxides have their largest values of T c when the 
Fermi surface is strongly influenced by Van Hove singu- 
larities. Then correlation effects can be estimated to be 
larger in the oxides and, finally, fullerides seem to have 
rotational disorder which would favour s-wave. Therefore, 
in principle, it could happen that in the oxides, going into 
the overdoped phase might lead to a crossover from d-wave 
to s-wave depending on the parameters. 



4 Conclusions 

In isotropic s-wave superconductors, the first nonadiabatic 
corrections to the e-ph interaction such as vertex and cross 
functions are strongly dependent on the momentum trans- 
fer q. In particular, small values of q leads to positive 
nonadiabatic corrections inducing an enhancement of the 
critical temperature T c Here, we have addressed the 
problem of the momentum dependence of the nonadia- 



batic corrections for a d-wave symmetry of the order pa- 
rameter. By introducing a model interaction in which the 
e-ph interaction is dominant at small values of q and the 
residual repulsion of electronic origin is instead impor- 
tant at larger momentum transfers, we have shown that 
also when the solution has d-wave symmetry, the inclu- 
sion of nonadiabatic corrections enhances T c compared to 
the case without corrections. Therefore in a strongly cor- 
related system, for which the e-ph interaction is mainly 
of forward scattering, d-wave superconductivity driven by 
phonons can be sustained by the nonadiabatic corrections 



A Analytical calculation of vertex and cross 
functions 

A.l Vertex function 



The evaluation of the vertex function given in the Eq.(33) 
follows basically the same lines and approximations made 
in Ref.[Q, the main difference being that here we refer to 
a two-dimensional system rather than a three-dimensional 
one. Making use of the linear model for the electronic dis- 
persion and considering the limit of T c /luq <C 1, we obtain 



l9(k-q)| ! 



£q-k+k' - 



+ U)Q — iuJ ri 



w 



0(e, 



q-k+k' 



lU! n 
-Eq-k+k') 



£q-k+k' + Wo 



Cq-k+k' 



LOO 



.(45) 



The main difficulty comes from e q -k+k' which, within our 
model, is 

e q -k+k' = v-p[q 2 + k' 2 + k 2 - 2qkcosa + 2qk' cos (3 

- 2kk'cos8} 1/2 -n, (46) 

where a, (3 and 8 are the angles between the directions of 
(q, k), (q, k') and (k, k') respectively. In the limit of small 
<7c the presence of (9-function in front of Py [Eq.(41)] and 
inside of the integral leads to |k| ~ |k'| and |k| ~ |q|. 
Therefore Eq. ( ff6| ) can be rewritten as follows: 

fq-k+k' — £q + v-pk-p [1 — cos a + cos /3 — cos 8] , (47) 

where we have taken |k| ~ kp. We can relate the angle (3 
to a and 8 by means of relation — 8 — a. Therefore the 
Eq. ([471) becomes: 



e q _k+k' — £q + wf^f [(1 — cos a)(l — cos ( 



sin a sin 8] . 
(48) 



If Q = [k - k'|/2fc F then cos0 = 1 - 2Q 2 and sinfl 

2Q\A ~Q 2 - 



e q-k+k' 



e q + EQ 2 (1- cos a) + EQ^1 - Q 2 sina . (49) 



Expanding cos a and sin a in powers of a and retaining 
only the lowest order term in a and Q we finally obtain: 



Cq-k+k' 



EQa 



(50) 
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For small q c we can replace 



9(q c - |k - q|) ~ 6(q c - k F y/2(l - cos a)) 

~9(2Q c -\a\). (51) 

At this point it is convenient to transform the momentum 
integration into an energy integration: 



The integrations over the energy and the angle are elemen- 
tary, the final expression of the cross function is however 
quite complicated: 



Pc(k,k') = U>oB(u> n , —U m ) 



1 



d 2 q AT r da f E ' 2 J 
—-L- = N / — / de 

(2tt) 2 °J-„2n J_ E/2 



(52) 



where we have used a constant DOS N(e) = No in the 
range —E/2 < e < E/2. The integration over the energy 
e and over the angle a can be performed analytically and 
the final expression of the vertex function in the limit of 
small Q c is given by 



r, <■>> , w ° 1 * ( 2 Q< E Q 

Pv(k, k ) = ujoB{uj ni oj m ) + ttt^-^t^ arctan 



2QcEy/l-Q 2 p{k,k f ) 

x {cos[r/(fc, k')]C(k, k') + sinfo(fc, k')]D(k, k')} 

A(u> n ,-u> m ) ~ B(w n ,-u> m )(u> n +uj m ) 2 ,_„s 
x i 1 i ' 

B are the 
— u> m . The function C, D, 77, p are 

1/4 



where the functions A and B are the same of Eqs. (|5_ 
and ([[5]) with u) m 
given by 



A(uj n ,u) m ) - S(w n) w OT )(w„ - u m ) 2 



T](k, k) = — — arctan 



where Q = |k — k'|/2fcp, Q c = q c /2kp, and 



(53) 



C(k,k') = b(k,k') 2 /arctan 



2|w n + o; m | 

£Q 2 

a+(fc, fe') 



(59) 
(60) 



A(uj n ,uj m ) = (uj n - u m ) 
+ arctan 



arctan | — — arctan — 
100 ) V 



arctan 



- arctan 



arctan 



"a_(fc,fc')~ 




+ arctan 


b(k,k') _ 



b(k,k') 
'2Q c -a + (k,k') 



b(k,k') 



2Q c -a_(fc,fc') 
b(k,k') 



(61) 



Lu + E/2 



LU Q + E/2 



(54) 



B(uj n ,Lu m ) = -(u) + E/2) 



A. 2 Cross function 



(lo + E/2) 2 + 2uj 2 n - uj n uj n 



D(k k>) _ 1 f [2Q c -a + (k,k')] 2 + b(k,k'r 
a-{k,k') 2 + b(k,k'f 



[(u + E/2)2 + u 2 



(55) 



a±(k, k') = 



a+(k,k') 2 + b{k,k') 2 

±p(k,k')cos[f](k,k')} 



b(k,k') = ^M= S in[ V (k,k')} 



The function Pc(k,k'), given by the Eq. (|34|) , can be ex- 
plicitly evaluated within the same scheme of calculation 
of the vertex function. In the limit of T c /loo <C 1 we have: 



A. 3 s- and e?-wave averages 



(62) 
(63) 
(64) 



Pc(k,k') 



l.9(k-q) 



wo 



In what follows we report the final expressions of the s- 
wave and d-wave averages of the vertex and cross functions 
defined in Eq.(f4l"|) for I — (s-wave) and I = 2 (<i-wave), 
where Py and Pq are given by Eqs. (|||) and (|5|), respec- 
tively. 



+ 



e q + lo - iuj n e q - uj q - iuj„ 

^(Cq-k-k') 9\— Cq-k-k') 



£q-k-k' + ^0 + £q-k-k' — ^0 + 



.(56) 



Py (u> n ,u) m ;q c ) 



A(bJ n ,U m ) - B(uj n ,U m )(u) n - LUn 



The electron energy e q _k-k' can be approximated for q c <C 
2kp as follows: 

e q -k-k' = Vf [q 2 + fc' 2 + k 2 — 2qk cos a — 2qk' cos /3 

1/2 _ 



w ° ^ / \ i t->/ arcsinQ c 
x tt^7 ^fiK,w m ,Q c ) +B(w„,w m ) — (65) 



P c (u; n ,uj m ;q c ) 

A(u n , -u) m ) - B(uj n , -u m ){u> n + u) m ) 2 



2kk'cos9} 1/2 - n 



k„ + uj m \ 



e q + E(l ~ Q 2 ) — - EQ^l - Q 2 a . (57) 



arcsin Q c 
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^rmvi^n, ^m] Qc) 



A(uj n ,ui m ) - B(uj n ,(j m )(u n - uj n 



u 
2EQI 



F 3 (LO n ,u m ,Q c ) + B{uj n ,uj m ){\ - 2Ql)^l-Ql 

(67) 



Pq 2 { UJ n, ^>m'i Qc) 

A(lJ n ,-U> m ) - B(LU ni -LU m )(LU n +uj m ) 2 



x ^k F ^n, u> m , Qc) + B(w n , -w m )(l - 2Q 2 c )^T~Ql 



2EQ 
where 

f 

Fl(u> n , U m , Qc) = I 

Jo 

n i L0 m , Qc. 



(68) 



dQ 



Qy/T^Ql 

r Q ° dQ 



: arctan 



2EQ C Q 



(69) 



1 



x {C(fc, k') cos[r/(fc, jfe')] - £)(jfe, k') sm[ V (k, k')]} (70) 
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jti / n , f Qc dQ ( 1 + 8Q 4 - 8 



x arctan 



2EQ C Q 



(71) 



T-i / n , [ Q *dQ/l + 8Q 4 

f4\U n , LO m , IJc) — J — I - 



1 



Q 2 ) P(k,k') 
x {C(k, jfe') cos[?7(fc, k')} - D(k, k') sin[7,(fc, k')}} (72) 



The functions C, D, 77 and p are precedently defined. 
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